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ABSTRACT 

Six different free fall missions that utilize the gravitational influence of a 
passing planet are presented. The missions are: minimum flight time for deep 
space probes, maximum inclination out of the ecliptic, maximum component of 
the velocity normal to the ecliptic, maximum distance out of the ecliptic within 
a specified time, minimum closest approach to the sun, and attaining a specific 
value of perihelion. Any mission objective that can be defined as a function of 
the heliocentric orbital parameters may also be treated. An iterative procedure 
is used to optimize the post encounter orbit for the desired mission. 

The patched conic assumption is used in that the actual flight path of a free 
fall vehicle is approximated by conic trajectories in finding the post encounter 
trajectory. 

Examples are included of missions utilizing the gravitational influence of 
the planet Jupiter. Without using Jupiter for an assist such missions would be 
difficult to perform with our present technology. 
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LIST OF SYMBOLS 


Notation 

r region of sphere of influence of target 

Vj asymptotic approach velocity vector of probe with respect to the 

sun as it enters r 

m mass of target planet 

/ip gravitational constant of planet = Gm 

M mass of sun 

u s gravitational constant of sun = GM 

| B | = BMAG miss distance 

I 

Vj asymptotic approach velocity vector of probe with respect to tar- 

get as it enters r 

— » 1 

V 2 asymptotic departing velocity vector of probe with respect to tar- 

get as it leaves r 

V 2 asymptotic departing velocity vector of probe with respect to sun 

as it leaves r 

V j heliocentric velocity of planet when probe enters r 

V 2 heliocentric velocity of planet when probe leaves t 

S radius of sphere of influence 

rP heliocentric position of planet when probe enters r 

‘ ‘ ♦ 

heliocentric position of planet when probe leaves r 
r p pericenter 

r a apocenter 
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GRAVITY ASSIST OPTIMIZATION TECHNIQUE 
APPLICABLE TO A 
VARIETY OF SPACE MISSIONS 


I. INTRODUCTION 

When a space vehicle is launched on a free fall trajectory in such a manner 
that it comes within the vicinity of a planet, the gravitational field of the planet 
can considerably alter the vehicle's orbit about the sun. The approach trajectory 
can be established in such a manner that the post encounter trajectory will allow 
the probe to complete many missions that otherwise could not be accomplished 
by a direct flight from Earth without extreme energy or payload penalties. For 
example, if it is desired to fly a solar probe mission, a direct flight must have 
sufficient energy not only to escape the Earth's gravitational field but also to 
negate the Earth's orbital velocity. If the probe is launched with a hyperbolic 
excess velocity of 11 km/sec, the closest it will come to the sun is .25 A.U. 1 
If, however, the gravitational influence of the planet Jupiter is used, it is then 
possible to impact the sun with approximately the same hyperbolic excess 
velocity. 

In a similar manner, many interesting missions may be accomplished which 
otherwise would be extremely difficult or expensive. It is the purpose of this 
report to develop a theoretical method of utilizing the gravitational influence of 
a planet in order to accomplish some of these missions. 


H. FLYBY-DETERMINATION OF POST ENCOUNTER CONDITIONS 

In determining the post encounter trajectory after it has passed through the 
gravitational field of a planet, it is assumed that certain pre-encounter condi- 
tions are known. 

Assuming a fixed transfer time from a park orbit about the earth to the 
planet or target we may obtain an approximate transfer trajectory to the target 
as in Reference 2*. This trajectory assumes both the earth and the target 
planet to be massless points. We obtain as a result the velocity vector (Vi) of 
the probe with respect to the sun at encounter with the target. 


♦Subroutine PLANET 
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We then take into consideration the gravitational field of the planet and 
utilize the patched conic assumption which states that at any given time the 
probe is under the influence of only one body. If m and M are the masses of the 
planet and sun respectively then the radius (S) of the sphere of influence is 
given by 


S 



R p 


where R p is the planet's distance from the sun. We have V t and knowing the 
injection date of the probe and its transfer time we can find V 1 the heliocentric 
velocity of the target. Reference to Figure 1 shows that V, = V, - V ? where 
Vj is the velocity of the probe with respect to the planet. 



Within the sphere of influence r, we assume that the probe's trajectory will be 
hyperbolic with respect to the planet (see Figure 2). We now make a basic and 
very critical assumption. We assume that Vi is now the velocity of the probe 
when it enters r and that we may choose the radius of closest approach of the 
probe without significantly altering the direction or magnitude of Vj. This is 
equivalent to assuming that while the direction and magnitude of Vi are fixed at 
the calculated value, we are still free to choose the point at which the vehicle 
pierces the surface of r (point A in Figure 2). This equivalence is illustrated 
in Figure 2: for a given asymtote and energy at entrance to r, the entire tra- 
jectory within r, including the RCA, is determined by the location of point A. 
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That our assumption is valid is indicated by Table 1 where = V x - V i lies 
along the S° vector (Figure 2). This table was compiled for both retrograde 
and posigrade flybys of the planet Jupiter and for a radius of closest approach 
varying from 100,000 to 10,000,000 km.* The values for RCA span a far larger 
region than we are interested in. 


Table 1 


Posigrade 
RCA x 10 5 km 

s° 

IvJ 

Right Ascension 
(DEG) 

Declination 

(DEG) 

(km/sec) 

i 

161.38295 

-.037895218 

13.863505 

5 

161.36287 

-.037826118 

13.873625 

10 

161.34038 

-.037764311 

13.882193 

50 

161.10998 

-.037458739 

13.910651 

100 

160.75975 

-.037179910 

13.922712 

Retrograde 




1 

161.43148 

-.037922481 

13.863560 

5 

161.48549 

-.037896001 

13.873828 

10 

161.53529 

-.037877017 

13.882264 

50 

161.82623 

-.037871506 

13.910774 

100 

162.11483 

-.037958746 

13.923062 


A convenient method for specifying the location of point A and the radius of 
closest approach is to select the magnitude of the miss vector (|b| = BMAG in 
Figure 2) and an angle p, measured between the B vector and a fixed vector in 
the planet's orbital plane. 

Now, the total effect of the encounter is simply to rotate Vj thru an angle y, 
as illustrated in Figures 1 and 2, since we must have 

Ivj = Ivj 


if energy is to be conserved. 


* Table 1 was generated using the refine option of the Quick Look Program (Reference 2). 
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|B I = BMAG 

\R 2 \ = S 


Figure 2 

We may readily calculate y if BMAG is specified: 

let a' = semi major axis of the hyperbolic trajectory in r. 

then 

M p S 

3 ‘ = 2/x p - S|v ;| 2 

from the energy equation. 

From Figure 2 and the properties of an hyperbola we have the eccentricity 

1 

e ~ cos e where e is the half angle between the asymptotes 
|§| _ BMAG 

and sin e = — — — - ; — (cos e) 

a e a 

BMAG 

Therefore tan e - ; — 
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Now 7 


77 


2e or 


7-77 



In order to compute V 2 we now need to define the second variable 0 which with 
BMAG will define the entrance point on r. We establish an orthogonal coordinate 
system consisting of the unit vectors S°, T°, and R° in the following manner: 

S° - — — lies along the incoming asymptote. Let k° be a unit vector perpen- 

I vj | 

dicular to the planet's orbital plane, then 

T° - k° x s°, T° lies in the plane of the planet's orbit 

R° = S° x f° 

» _ < ^ 

Figure 3 shows V x , V 2 , and B in relation to this coordinate system. 



Figure 3 

Since B _L S° , B lies in the plane formed by T° and R°. Also, since the hyper- 
bolic orbit lies in a plane B , and V 2 must lie in the same plane. We define 
our second variable \jj to be the angle between T° and B measured from T° to 
R°. Now with reference to Figure 4 we may decompose V 2 into its components 
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along S\ T°, R° and obtain 

V j - !vj(s° cos y ~ sin y [r° sin \p + T° cos v//]J . 

For a first approximation we assume_the time spent in r to be negligible and 
therefore we let Vj r vJ,thenV 2 = vj + (the heliocentric departure 
velocity). Similarly if R 2 is the position vector of the departing probe we may 
assume R 2 = where R ^ is the heliocentric position of the planet correspond- 
ing to V ’j . 

We now use the approximate values R 2 and V 2 of the probe to obtain a more 
accurate estimate of the post encounter position and velocity. As before we have 

1 

Co S € ~ ~ 

e 


and Cos 6 




+ Cos 2 e 


and with reference to Figure 2 we see that 


V • V 

V 1 V 2 


| V 9 | Cos ( tt - 2 e) 


- iVj I V 2 I Cos 2 e 


then e 


\ 1 + Co 


s 2 e 


6 





2 


e 




p - a ( 1 - e‘ 2 ) - the semi latus rectum. To find the true anomaly at exit from 
r we use the radius equation and obtain 


e 


Cos - 1 



S 


The time from periapsis to true anomaly Q is calculated,* which is one half the 
total time spent in r. With this time, we may step forward and backwards along 
the heliocentric orbit to obtain new values for V j, V 2 ,R j, andR 2 : We then 
recalculate a 1 , 7, S°, T°, R°, V 2 , p and 8 exactly as we did before. We re- 
peat this iterative procedure until time from periapsis is within a predetermined 
tolerance. Then V 2 = V 2 + V 2 using the values from the last iteration only. 

We now need R 2 , the heliocentric position of the probe when it leaves r. To 
calculate this we first calculate a unit vector B° along B. From Figure 4 


B° = Cos 0 T° + sin 0 R° 
and from Figure 2 we find 

R 2 = S [b° Cos ( 8 + tj/2 ~ e) + S° Cos ( 6 - e)] 
then r 2 = r 2 + r 2 

From R 2 and V 2 we may now calculate the post encounter heliocentric 
trajectory. From the area integral 

r 2 x v 2 = h 2 


* Reference 2, Subroutine TCONIC . 
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1 _ 

The Hamilton Integral for the post encounter trajectory is V 2 = — H 2 * ( R 2 + e ) 

P 

where e is a vector in the direction of perigee. 

Upon cross multiplying on the right by H 2 we obtain 

1 - 

e = V 2 x H, - R°, 

PS 2 2 


then e 


also 


H 2 • H 2 

Ps 


(1 - e 2 ) 


and from this we obtain a. The true anomaly is easily obtained from the radius 
a (l - e 2 ) 

equation r — L — and the inclination by taking the dot product of H, 

1 + e Cos 9 1 

and a unit vector perpendicular to the ecliptic. We now have the standard 
elements needed to compute the total trajectory. 

HI. POST ENCOUNTER 

Maximum distance out of the ecliptic within a given time 

Appendix A describes a procedure for obtaining the maximum distance out 
of the ecliptic that any given post encounter trajectory attains. When the period 
of the orbit is large, it may not be realistic for a mission to be based on reach- 
ing such a distance. A more useful mission might be to attain the maximum 
distance possible within a specified time (At) from post encounter. 

From At and the position and velocity vectors (R and V) at departure 
from t, the position and velocity (R t , V t ) of the probe at time At can be 
computed.* Figure 5 shows the orbit with respect to the line of intersection 
with the ecliptic plane where the points P and P' represent the maximum 
distances out of the ecliptic as discussed in Appendix A. 

By defining u° to be a unit vector along the line of intersection it is easily 
seen that 


Reference 2, Subroutine STEPD. 



D 


1 

2 



D i = 
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If i is the inclination of the orbit to the ecliptic then the distances out of the 
ecliptic are 


d - D sin i 
dj = Dj sin i 

we now have computed four distances out of the ecliptic: the two maxima 
associated with p and p', and d and di. To determine which of these distances 
is the maximum reached within A t we define k° to be a unit vector perpendi- 
cular to the ecliptic and form the dot products V • k° and V t * k° . If these 
quantities are both positive or negative then the maximum distance is the 
larger of d , dj. If the dot products are of opposite sign then the proper solu- 
tion is one of the distances associated with P and P'. The sign of V • k° 
determines the correct one. Finally if either dot product is zero then we are 
already at the maximum. 

The distance thus found is the functional value to be used in the optimization 
process. 
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Minimum Time for Deep Space Probes 


In this option we compute the time from post encounter to a given heliocen- 
tric distance (x A.U.’s). Using the post encounter elements, apogee is calculated 

by 

r a = a (l + e) 


For this option we therefore must have 

r > x A.U. 

a — 

Knowing the true anomaly at departure enables us to calculate* the time t Q 
from periapsis to that value of the anomaly. The true anomaly of the probe at 
r - x A.U. is 

a (l - e 2 ) - r 
6 - Cos -1 Fe 

and from this we may calculate the time (ti) from periapsis to x A.U. The 
functional value we then wish to minimize is 


Maxi mum Inclination to the ecliptic 

From the position and velocity vectors (R and V) of the post encounter 
trajectory we obtain the angular momentum vector H = R * V . Ifk°isa unit 
vector perpendicular to the ecliptic then 

H * k° 

= Cos i where i is the inclination 

I H | 

The maximum value of i occurs at 90° or Cos i = 0* The functional value for 
the optimization process is |Cos i|. 

Maxim um velocity component normal to the ecliptic 

For this option we easily obtain that component of the velocity vector which 
is normal to the ecliptic from V • k° , where V and k°are defined as above. In 


* Reference 2, Subroutine TCONIC. 
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the ideal case this would occur at perigee and therefore the post encounter 
trajectory should have a true anomaly near zero. 

Minimum perihelion distance 


Perihelion (r p ) is determined from r p = a (1 -e). For this option it is 
desirable for the probe to be moving toward perigee so that perihelion is 
reached within a reasonable time. The optimization routine therefore constrains 
the true anomaly to be greater than 180 degrees. 

Specific perihelion distance 

In this option it is desired to determine a post encounter trajectory which 
makes r p equal to a specified value. Again true anomaly is constrained to be 
greater than 180 degrees. The functional value to be minimized is the absolute 
value of the difference between the desired and actual values of r p . 

Additional Missions 


The above options were analysed as examples of the class of mission ob- 
jectives that are useful in a gravity assist program. 

Further options that can be defined as functions of the post encounter 
trajectory parameters may be included and treated in a manner similar to the 
above missions. 

IV. OPTIMIZATION PROCEDURE 

As previously described, BMAG and \p completely determine the post en- 
counter trajectory for a given Vi. This orbit is then used to obtain one of the 
post encounter missions described above. 

However, for any given mission, the original independent variables wall not 
be those which give the optimum for the mission desired. To accomplish the 
optimization of the particular mission, a multi-variable nonlinear iterator* is 
used w'hich systematically alters the independent variables, calculates a new 
post encounter trajectory and a new functional value for the iterator to use in 
the next iteration. The functional value is that quantity to be maximized or 
minimized depending on the mission chosen. The iterating process continues 
until an optimum value is achieved. 


* Reference 3. 
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V. EXAMPLES 


In Tables 2 through 7, examples of each mission is given for five and six 
hundred day transfer times to the planet Jupiter. It is seen that some trajectories 
may be used for more than one purpose. For example the 500 day transfer for 
attaining a perihelion of .25 A.U. also may be used for out of the ecliptic studies 
since it reaches a distance of .9 A.U. from the ecliptic. As one might expect, 
missions out of the ecliptic show a certain symmetry in that the probe can pass 
above or below the planet with one trajectory being essentially an inversion of 
the other. 

Tables 8 and 9 show that for a specific perihelion and for maximum inclina- 
tion out of the ecliptic there is a range over which p may vary and still permit 
a successful mission to be performed. In Table 8, for between ±15° the peri- 
helion distance of .25 A.U. is reached by varying BMAG. Table 9 shows that 
high inclination is attained along with useful distances out of the ecliptic for the 
ranges 


10° < \p<_ 60° 

- 60° < \p < - 10° 

For \p - ±10°, the probe will reach only about l/2 A.U. out of the ecliptic before 
impacting the sun. However, for a multipurpose mission this may be acceptable. 
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Maximum Inclination to the Ecliptic 
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Table 4 
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Table 5 
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Table 6 
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Desired Perihelion Distance of .25 A.U. 
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APPENDIX A 


Maximum Distances out of the Ecliptic 


When the post encounter trajectory takes the probe out of the ecliptic plane 
it would be meaningful for certain missions to know the maximum distance out 
of the ecliptic that the probe attains. Clearly this implies that the trajectory is 
an ellipse. Such a trajectory with a relatively small period would be extremely 
useful in studying those regions of space above and below the ecliptic. 

Figure 6 shows a typical post encounter trajectory which is inclined to the 
ecliptic by an angle i . The vector e is directed toward perigee; H = R x V 
and is perpendicular to the orbital plane and k° is a unit vector perpendicular 
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to the ecliptic. 9 e is the angle formed by the vector e and the vector k° x H 
which lies along the line of intersection of the two planes. The angles i and 9 
are easily found from the post encounter trajectory as follows: 


Cos i 



Cos 9 

e 


(k° x H) • e 

12° x H| | e| 


To compute the maximum distances both above and below the ecliptic it is 
convenient to consider a coordinate system with origin at the center of the ellipse 
and with the X axis along e . Referring to Figure 7 we compute the maximum 
distance from the ecliptic at the point (X o , Y 0 ). 



Figure 7 
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From the post encounter trajectory, a and e - |e| are already known. The 
equation of an ellipse is 


a \* 


±b 1 


where 


and 


/r 7 " 


dy 

d x 


at (X o , Y o ) we must have 


b 


Xb 


r.l \ ~' A 


1 " 


X — 1 - 


X 2 

o 




± Tan 9 


so that 


then Y 


a 4 Tan 2 9 


b 2 + a 2 T an 2 9 e 
a 2 Tan 2 9 


X 


1 - e 2 + Tan 2 9 


(i 


( X 2 

) 


substituting for b 


where we consider only the absolute values of X Q and Y o . From Figure 7 

Y 


D' 


Cos 9 


Sin 9 

e 


x 


o 


IT 

g + a e 


X > 0 

o — 
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thus D 


X Sin 0 - g Sin 8 + ae Sin 8 

o e e e 

Y c 

also Tan 90-0 = cot 8 - — ; Y > 0 

e e f 9 o — 

Sin 8 ' 

then f = Cos 8 

e 

Y o Sin 2 8 e Y o 

finally D = D‘ + D" = X Sin 8 7 + ae Sin 8 + 7 ; 

0 e Cos t? e Cos 6 


or D - X Sin 0 + Y Cos 6 + ae Sin 9 (1) 

o e o e e ' ' 


for the distance on the opposite side of the ecliptic we obtain similarly 


D, = Y Cos 8 + X Sin 8 - ae Sin 8 (2) 

1 o e o e e 

Finally, the distances out of the ecliptic are 


dj - D Sin i 


d 2 - D 2 Sin i . 

The distance computed using equation (1) will be the maximum distance 
(d 1 ) attained by the probe. That distance obtained from equation (2) is the 
maximum distance on the other side of the orbit (d 2 ). If k° • e > 0 then d t 
is positive and above the ecliptic while d 2 is below the ecliptic; the opposite is 
true if k° • e is less than zero. 


25 



APPENDIX B 


BLOCK DIAGRAM OF COMPUTER PROGRAM 


QUICK LOOK PROGRAM calculates V| and 
time from 1950 


I 


SUBROUTINE FXMISS 







calls FLYBY to calculate 

calculates position and velocityof target 



the post encounter trajectory 

planet at given time 



calls FUNCTN to calculate 

calls FCTS for initial printout 

— ► 

to 

u 

the function to be optimized. 


LL 




LU 

Z 


initialize for the iterator 


h- 

3 

O 

OH 

CO 

calculates the true anomaly 
constraint if true anomaly 



3 

cO 

is to be negative (min. 
perhelion) 

calls MINMAX (iterator) 




calls FCTS for final printout 

1 
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